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1.  INTRODUCTION 


The  spin  of  all  projectiles  decreases  along  the  flight  path  because 
of  the  aerodynamic  torque  arising  from  the  skin  friction  on  the  exterior 
surface  of  the  casing.  The  spin  of  a liquid-filled  shell  decreases  more 
rapidly  than  one  with  a solid  payload  because  angular  momentum  is  trans- 
ferred from  the  projectile  casing  to  the  liquid  payload  while  the  liquid 
is  being  spun  up.  I.iquid  spin-up  is  produced  by  tlie  skin  I'riction  at 
the  casing- liquid  interface.  We  define  the  time  in  which  the  liquid, 
starting  from  zero  spin,  achieves  a state  of  substantially  solid-botly 
rotation  to  be  the  spin -up  period. 

A study  of  the  liquid  motion  during  spin-up  is  imjiortant  for  several 
reasons.  First,  the  frequencies  of  free  oscillation  of  the  liquid  change 
during  spin-up.  These  frequencies  are  needed  to  analyze  the  flight  sta- 
bility of  a 1 Iquid- fi 1 led  shell.  Second,  spin  decay  can  decrease  the 
gyroscojiic  stal^ility  factor  of  a liquid-filled  shell  to  an  unacceptable 
level.  Third,  the  projectile  spin-decay  jiroccss  may  ]X)ssibly  lead  to 
fluid  dynamic  instabilities  in  the  liquid  payioad  which  could  affect  the 
shell  motion.  In  rhis  report  oniy  the  spin-decay  process  is  studied, 
i'he  effect  of  s]un  decay  on  the  frequencies  will  be  reported  sejiarately. 

The  basis  of  tl\e  present  work  is  a tiicory  dcvelo])od  by  Wedemeyer' 
which  describes  spin-up  from  rest  of  a liquid  in  a fully-filled  cylin- 
drical cavity,  iiu  s theory  accounts  for  a secondary  flow,  formed  in  the 
c:ivity  as  a result  of  tlic  endwall  boundary  layers,  which  controls  the 
spin-up  process.  Wedemeyer  developed  an  equation  for  the  azimuthal 
velocity  during  sjun-up  and  solved  for  the  velocity  in  closed  form  by 
neglecting  viscous  diffusion  terms,  lie  used  this  solution  to  derive 
approximate  expressions  for  the  rate  of  change  of  angular  momentum  of 
the  liquid  and  calculate  spin-decay  of  liquid-filled  shell. 

In  this  report  wc  also  calculate  spin-decay,  but  our  analysis  dif- 
fers fi'om  Wedemeyer' s in  several  respects.  Hirst,  we  adopt  Wedemeyer' s 
azimuthal  velocity  equation  describing  the  spin-up  process,  but  retain 
the  viscous  diffusion  terms  he  neglected.  Wc  determine  the  fluid  angu- 
lar momentum  numerically  from  finite-difference  solutions.  Second,  our 
analysis  includes  a teclmique  for  treating  the  non- impuls i ve  spin-up 
process  taking  place  in-bore  the  gun.  This  is  used  to  determine  the 
initial  conditions  for  the  subsequent  spin-decay  calculation  in  free 
flight.  The  accuracy  of  this  technique  is  evaluated  by  comparing  its 
results  with  those  of  finite-difference  solutions  of  ihe  Nav i er-Stokes 
equations;  it  is  Judged  to  be  sufficient.  Third,  tlie  aerodynamic  spin 
damjiing  effect  is  included  in  our  analysis  so  the  projectile  spin  can 
be  accurately  predicted  throughout  the  entire  flight. 


1,  E.  //.  Wadenieijev,  "The  Uns beady  Flow  Within  a Spinning  Cylinder," 

J . Fluid  Me  oh,,  Vol.  "0,  Part  lOC'l,  pp.  'S83-399,  Also  see  BHL 
Report  Aberdeen  Proviruj  Ground,  MD,  Oetoher  1963.  (AD  '131896) 
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Spin-decay  predictions  obtained  with  these  techniques  are  compared 
with  measurements  taken  in  test  firings  of  the  XM687  shell  using  a solar 
aspect  sensor  and  telemetry  technique^.  The  predicted  spin  decay  is  in 
excellent  agreement  with  the  measurements  except  for  certain  cases  which 
show  anomalous  spin  behavior;  these  cases  are  discussed  in  a separate 
report^ . 


II.  PROJECTILE  SPIN-DECAY  EQUATIONS 

Consider  a projectile  containing  a liquid-filled  cylindrical  cav- 
ity. Translational  motion  of  the  projectile  in  the  gun  imparts  spin  to 
it  because  of  the  barrel  rifling.  Spin-up  of  the  liquid  begins  in-bore 
and  continues  after  the  projectile  exits  the  gun  tube.  In  free  flight 
the  spin  rate  of  the  casing,  p,  begins  to  decrease  from  its  launch 
value,  p^,  because  of  moments  produced  by  two  shear  forces:  1)  spin- 

decelerating  moment  due  to  air  shear,  and  2)  liquid  shear  moment, 

Mbiq*  which  causes  angular  momentum  to  be  transferred  between  the  casing 

and  the  liquid  payload.  The  projectile  spin  rate  is  determined  from  the 
equation  of  motion 


I dp/dt  = M,  . + M.  ; 

2 Liq  Aero  ' 

(1) 

where 

P V2  „ 

^ero=  2 ' 

(2) 

P 

^Liq  * • 

(3) 

The  quantity  ^2  axial  moment  of 

inertia  of  the  empty  projectile. 

V is  the  projectile  velocity,  S is  the  maximum  projectile  cross-sectional 
area,  I is  the  projectile  diameter,  p is  the  air  density,  C„  is  the 

p 

projectile  roll  damping  coefficient,  and  L is  the  liquid  angular  momen- 
tum. Since  is  negative  over  the  whole  trajectory  it  causes  p to 


2.  A,  Mark  and  W.  //.  Mermagen,  "Measurement  of  Spin  Decay  and  Insta- 
bility of  Liquid-Filled  Proyeatiles  Via  Telemetry,"  BRL  Memorandum 
Report  2333,  Aberdeen  Proving  Ground,  MD,  October  1973.  (AD  771919) 

3.  C.  W.  Kitchens,  Jr.,  and  R,  Sedney,  "Conjecture  for  Anomalous  Spin 
Decay  of  the  155mm  Binary  Shell  (XM697),"  BRL  Memorandum  Report  (in 
preparation) . 
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decrease.  ^L^q>  other  hand,  is  negative  at  launch  and  in  the 

early  portion  of  the  trajectory,  but  positive  late  in  the  flight.  The 
liquid  spin-up  process,  with  dL/dt  > 0,  takes  place  during  the  in-bore 
and  early  free-flight  phases.  The  moments  acting  on  the  casing  during 
these  phases  are  shown  in  Figures  la,  b.  While  the  casing  spin  rate 
decreases,  liquid  spin-up  proceeds  until  the  liquid  payload  approaches  a 
state  of  .solid-body  rotation.  Solid-body  rotation  is  never  achieved 
because  continues  to  decrease  the  projectile  spin  rate.  The 

liquid  will  eventually  be  "overspun"  relative  to  the  instantaneous  pro- 
jectile spin  rate;  that  is,  the  local  liquid  angular  velocity,  w = v/r, 
will  be  greater  than  p,  where  v is  the  liquid  azimuthal  velocity  and  r 
is  the  radial  coordinate.  The  liquid  moment  then  reverses  direction; 
i.e.,‘ dL/dt  < 0,  thus  opposing  the  direction  of  the  aerodynamic  moment. 
This  is  illustrated  in  Figure  Ic.  During  this  phase  the  liquid  acts  to 
oppose  further  spin  decay  of  the  projectile. 


We  must  know  M.  and  M,  . as  functions  of  p and  t in  order  to  be 
Aero  Liq  ' 

able  to  predict  the  spin  rate  from  Eq.  (1).  The  first  is  obtained  from 
flight  measurements  using  a projectile  with  a solid  payload;  M . is 

Ij  1 t| 


determined  by  solving  the  equations  of  motion  for  the  liquid.  The 
boundary  conditions  for  these  equations  involve  the  projectile  spin  rate 
which  is  the  quantity  we  wish  to  determine.  If  the  fluid  motion  is 
known,  M^iq  found  by  integrating  the  shear  stress  acting  on  the 

sidewalls  and  endwalls  of  the  cavity  to  obtain  the  resultant  torque.  An 
equivalent  procedure,  which  is  more  convenient  for  the  present  purpose, 
uses  Eq.  (3).  I,(t)  is  defined  by 


h(t)  = p/J/r^v  (r,z,t)  dedrdz 


(4) 


where  the  integral  is  taken  over  the  cavity  volume;  r,  0,  z are  non- 
rotating cylindrical  coordinates  with  z along  the  axis  of  the  projec- 
tile, and  p is  the  liquid  density. 


In  our  analysis  we  neglect  the  projectile  yawing  motion  and  a.ssumo 
that  the  spin-up  flow  is  axisymmetric.  Wedemeyer's  spin-up  theory'  is 
used  to  determine  v. 


III.  CALCULATION  OF  LIQUID  ANGULAR  MOMENTUM 
A.  The  Wedemeyer  Spin-Up  Model 

Wedemeyer'  considered  the  problem  of  a circular  cylinder  filled 
with  liquid  which  is  initially  at  rest.  At  t = 0 it  is  given  an  impul- 
sive angular  velocity  about  its  longitudinal  axis  and  maintained  at  that 
value.  He  sought  to  determine  the  unsteady  flow  for  all  t,  until  the 
spin-up  process  ends  with  the  fluid  rotating  uniformly  as  a rigid  body. 
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The  flow  region  is  divided  into  two  parts:  (1)  the  thin  I-kman  boundary 

layers  on  the  cylinder  end  walls;  and  (2)  the  remainder  of  the  flow, 
called  the  core  flow.  He  finds  that  the  spin-up  process  is  dominated  by 
the  Ekman  layers.  They  create  centrifugal  pumping  which  causes  secondary 
flow  in  the  core.  This  flow  enters  the  Ekman  layers  ^.t  r(0  < r < r*), 
is  spun  up,  and  then  ejected  back  into  the  core  at  r > r*,  where  r*  is 
the  radial  position  of  a shear  layer  propagating  inward  from  the  side- 
wall.  This  spinning  fluid  is  then  carried  back  into  the  interior  of 
the  cylinder.  This  mechanism  of  imparting  rotation  to  the  fluid  is  much 
more  efficient  than  viscous  diffusion  alone. 

Wedemeyer  developed  an  approximate  equation  for  the  core  flow 
through  an  order  of  magnitude  ar.alysis  of  the  axisymmetric  Navier-Stokes 
equations.  In  the  core  v(r,2,t)  = v^(r,t).  The  tangential  momentum 

equation  describing  the  core  flow  becomes 

3v  /at  + u (avyar  + v /r)  = v[a^v  /ar^  + a(vyr)/ar]  ; (.S) 

Cl  O 0 O C O 

where  u^  and  v^  are  the  radial  and  tangential  velocity  components  in  the 

core  flow  and  v is  the  liquid  kinematic  viscosity.  The  analysis  shows 
that  u^  and  the  axial  velocity  w^,  the  secondary  flow  referred  to  above, 

are  small  compared  to  v^,  but  not  negligible,  and  that  u^  and  v^  are 

independent  of  z.  The  pressure  P is  given  by 

c 

3P^/ar  = pv^2/r  . (6) 

Conservation  of  mass  relates  u to  the  outward  radial  mass  flow  in  the 

c 

Ekman  layers.  An  assumption  is  then  made  for  the  relation  between  the 
Ekman  layer  mass  flow  and  v by  interpolating  between  the  known  results 
for  t -►  0 and  t ®.  ^ 

Wedemeyer  obtains  two  relationships  for  u^  depending  on  whether  the 
endwall  boundary  layers  are  expected  to  be  laminar  or  turbulent: 

u^  = -0.443  (a/cjRe"'^  (rp-v^)  , for  Re  < 3x10^  (7) 

and 

u^  = -0.035  Ca/c)Re'^^'^  (rp-v^)®'^^/(ap)^'^^  , for  Re  > 3x10^  ; (8) 

where  c is  the  cavity  half-height,  a is  the  cavity  radius,  and  the 
Reynolds  number  is 

Re  = p a^/v  . (9) 
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In  his  formulation  of  the  spin-up  problem  he  assumes  p is  constant,  but 
in  his  application  of  this  to  spin  decay  p is  a function  of  time.  We 
shall  also  take  p to  be  a function  of  time,  both  in  free  flight  and  in- 
bore the  gun.  A discussion  and  justif i'"’tion  of  the  latter  situation  is 
given  in  Section  III.C. 


The  initial  and  boundary  conditions  for  Eq.  (5)  are 
v^(r,t)  = 0 for  t < 0, 

V (0,t)  = 0 and  v,(a,t)  = ap  for  t > 0 

C (m- 


(loj 


During  projectile  in-bore  travel  p(t)  is  specified  using  the  spin  his- 
tory calculated  from  an  interior  ballistics  trajectory.  In  free  flight 
p(t)  cannot  be  independently  specified;  it  must  be  calculated  from  Eq. 
(1)  including  both  the  aerodynamic  and  liquid  moments. 

The  total  angular  momentum  of  the  liquid  within  the  cylinder  can  be 
expressed  using  Eq.  (4)  with  v equal  to  v^,  yielding 

L = 47TCP  /^r^v^dr  . (11) 

B.  Computation  Procedures 


The  authors  have  employed  two  procedures  using  Eq.  (1)  to  calculate 
spin-decay  of  liquid-filled  shell.  Method  I is  based  on  a numerical 
solution  of  Eq.  (5);  Method  II  is  based  on  Wedemeyer's  approximate 
expressions  for  rate  of  change  of  angular  momentum.  A third  method, 
Method  III,  uses  a Navier-Stokes  calculation  to  predict  in-borc  liquid 
spin-up. 

In  Method  I calculations  are  performed  by  simultaneously  solving 
Eqs.  (1),  (5)  and  (11),  using  (2)  and  (3).  Eq.  (5)  is  solved  by  a 
second-order  accurate  finite-difference  technique  described  by  Sedney 
and  Gerber**,  In  Method  1 we  rewrite  Eq,  (1)  as 

I dp/dt  = -dL/dt  + p f(t)  ; (12) 

2i 


where 


f(t)  = p^V  S n (13) 

P 

is  a known  function  of  time  depending  on  the  projectile  shape,  veloc- 
ity and  the  trajectory.  We  assume  that  all  variables  are  known  at  time 


4,  R.  Sedney  and  N.  Gerber,  '^Visaoua  Effects  in  the  Wedemeyer  Model  of 
Spin-Up  from  Rest,"  BRL  Report  (in  preparation) . 
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either  through  initial  conditions  or  subsequent  calculations  and  we 
wish  to  determine  p and  L at  time  We  can  express  hq.  (12)  by  a 

second-order  accurate  finite- difference  representation 

I (P  .,-P  = -(L  ,-L  )/At  + [p  ,f(t  J + p f(t  )]/2  (14) 

z^*^n+l  ^ n+1  n-^  *^*^n+l  n+1  ‘n  ^ 

and  solve  for  p , , to  obtain 
*^n+ 1 


P„(I,  • it  f(t„)/21 


L ♦ L 
n+1  n 


n+1 


[I^  - At  f(t„^p/2] 


(15) 


hq.  (15)  relates  the  spin  and  angular  momentum  at  t^^j  with 
f(tn^j),  Pj^,  At  and  1^  prescribed. 


The  fluid  angular  momentum  must  be  calculated  from  liqs.  (5) 

and  (11)  using  the  sidewall  boundary  condition  v^(a,t^^^j)  = ap(t^^^j). 

These  equations  give  us  a second  relationship  between  and 

which  is  very  complicated.  This  is  combined  with  Eq.  (15)  to  provide 
two  equations  for  P^^^^  ^n+1  solved  numerically.  We  adopt 

an  iteration  procedure  on  to  retain  second-order  accuracy  in  the 

solution.  The  procedure  employed  is  to  guess  using  linear  extrapo- 

lation based  on  values  at  previous  times  and  solve  Eq,  (!5)  for  Pj^^j- 
Eq.  (5)  is  then  solved  using  p^^j  as  a boundary  condition  to  yield 
V (r.t  ,)  The  angular  momentum  , is  calculated  and  compared  with 
the  guessed  value  If  they  differ  by  more  than  a set  tolerance,  t , 

the  calculation  is  repeated  using  as  the  guess  in  Eq.  (15).  I’his 

iteration  process  is  repeated  until  the  desired  tolerance  is  achieved. 

A value  of  e = L xlO“^  was  used  in  all  calculations  described  in  this 
o 

report  and  convergence  was  achieved  in  one  to  three  iterations.  This 
procedure  can  treat  either  laminar  or  turbulent  endwall  boundary  layers 
through  the  specification  of  Eq.  (7)  or  Eq.  (8)  for  in  Eq.  (5). 

Method  II  is  a simpler  procedure,  based  on  an  ordinary  differential 
equation  for  angular  momentum  which  is  solved  simultaneously  with  Eq. 
(12).  Equations  for  both  laminar  and  turbulent  cases  were  derived  by 
Wedemeyer^,  assuming  that  Eq.  (5)  is  valid  for  non-constant  spin  rate. 

In  the  laminar  case,  Eqs.  (5)  and  (7)  are  combined  and  integrated  in  r 
to  form  an  equation  for  the  time  rate  of  change  of  liquid  angular  momen- 
tum. This  equation  contains  a viscous  term  representing  the  moment  due 
to  the  shear  stress  along  the  sidewall.  Wedemeyer  approximated  the 
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shear  stress  using  an  analytical  "inviscid"  solution  obtained  by  neglect- 
ing the  viscous  terms  on  the  right-hand  side  of  Eq.  (5) . This  approxima- 
t ion  is  not  valid  late  in  free  flight  when  the  liquid  is  overspun  rela- 
tive to  Lhe  shell,  because  it  does  not  account  for  the  reversal  of  the 
liquid  moment.  Wedemeyer's  equation  for  angular  momentum  in  the  laminar 
case  is 


dL/dt  = - ,886  (a/c)(Re^)‘^(p/p^)*^  (Lp^'V^ 


8(Re^)-^  PLq  [1  - 


(16) 


where 


L = irpcp  a** 
o 


and 


Re  = p a^/v 
o ‘^o 


(17) 


L is  the  reference  angular  momentum  for  a fluid  in  solid-body  rotation 
o 

with  V - p r,  p is  the  instantaneous  shell  spin  rate  when  it  exits  the 
gun  muzzle  and  Re^  is  the  launch  Reynolds  number. 

For  the  turbulent  case,  Wederaeyer  used  the  same  procedure  with  Eqs. 
(5)  and  (8)  and  also  approximated  the  sidewall  shear  stress  using  the 
"inviscid"  solution  for  the  laminar  case.  An  additional  assumption  was 
required  to  evaluate  an  integral  which  depended  on  the  shape  of  the 
velocity  profile;  he  assumed  the  velocity  profile  could  be  approximated 
by  the  inviscid,  laminar  solution.  Under  these  assumptions  the  expres- 
sion for  the  turbulent  case  becomes 

dL/dt  - 0.035  (a/c)(Re^)"^/^  (p^L^) (p/p^)^/^  [ 1-Lp^/ (L^p) ] 


r.  /I  .1-13/5  .r^Po^^oP  „8/5  ,,  ..-3/10 
^ [LPq/LqP)]  4/q  X (1-X)  dX 


- 8CRe^)-^  P^o 

The  terms  on  the  right-hand  side  of  Eqs.  (16)  and  (18)  represent  the 
moments  acting  on  the  fluid  due  to  the  shear  stress  at  the  sidewall  and 

endwalls.  The  first  group  of  terms  in  each  equation,  containing  Re^  ^ 

and  Re^”^^^,  respectively,  represent  the  moment  due  to  the  endwalls. 

The  terms  containing  ^ represent  the  moment  due  to  the  sidewall. 

The  sum  of  these  two  moments  equals  the  time  rate  of  change  of  angular 
momentum  of  the  liquid. 
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To  calculate  the  shell  spin  decay  in  the  second  procedure,  Eq.  (16) 
or  Eq.  (18)  is  solved  simultaneously  with  Eq.  (12).  Although  simpler 
than  Method  I,  Method  II  has  the  drawback  that  the  approximations 
involved  in  evaluating  the  shear  stress  at  the  sidewall  (and  describing 
v^^(r,t)  in  the  turbulent  case)  are  not  accurate  during  the  late  free- 

flight  phase.  Fig.  Ic,  where  dL/dt  < 0.  We  have,  however,  used  both 
procedures  and  compared  their  predictions  against  the  measured  spin 
decay  of  liquid-filled  XM687  shell. 


C . In-Bore  Spin-Up  Effects 


The  first  step  in  calculating  projectile  spin  decay  is  to  determine 
the  amount  of  angular  momentum  achieved  by  the  liquid  during  the  in-bore 
spin-up  process.  This  prescribes  the  initial  conditions  for  the  subse- 
quent spin-decay  calculation.  Karpov^  carried  out  experiments  with  20mm 
M56  shell  to  study  the  spin-decay  process  and  instability  during  spin- 
up.  He  presented  spin  histories  for  firings  in  a vacuum  by  correcting 
for  spin  deceleration  due  to  the  aerodynamic  moment.  He  found  that 
higher  viscosity  liquids  reach  "constant  spin"  sooner  and  at  a level 
considerably  different  from  that  predicted  by  conservation  of  angular 
momentum.  He  thus  was  able  to  infer  the  percentage  of  the  reference 
angular  momentum  that  was  acquired  while  the  projectile  was  still  in 

the  gun.  L^,  Eq.  (17),  the  maximum  angular  momentum  that  can  be  achieved 


by  the  liquid  is  never  attained  in  actual  firings  because  of  the  short 

time  in  the  barrel  and  the  spin-decay  of  the  casing.  For  a 70"b-filled 

case  with  Re  = 6500  he  found  that  24%  of  L was  achieved  in-bore, 
o o 

For  higher  Reynolds  number  cases  the  percentage  was  smaller;  e.g., 

with  Re  = 6.5x10®  the  angular  momentum  at  muzzle  exit  was  10%  of  L . 
o " o 

While  these  actual  percentages  will  not  apply  to  other  shell  or  gun 

tubes,  they  illustrate  that  significant  liquid  spin-up  can  occur  in-bore 

the  gun. 


Wedemeyer^  neglected  the  in-bore  spin-up  process  in  predicting 
spin  decay;  the  angular  momentum  calculations  began  as  the  shell  exited 
the  gun.  Also,  Wedemeyer's  equation,  Eq.  (5),  is  based  on  the  assumption 
of  impulsive  spin-up  from  rest.  The  spin-up  process  actually  begins 
in-bore  from  a non-impulsive  start.  As  the  shell  accelerates  the  rota- 
tion rate  increases,  reaching  p„  at  muzzle  exit.  By  the  time  the  shell 

reaches  the  muzzle  the  liquid  may  have  a significant  portion  of  L^, 
especially  for  smaller  Re^.  Since  part  of  the  liquid  is  already  rotat- 
ing by  the  time  that  the  projectile  reaches  the  muzzle  exit,  the  spin 


5.  B.  G.  Karpov,  "Dymmiae  of  Liquid-Filled  Shell:  Inetability  During 

Spin-Up,'*  BRL  Memorandum  Report  1629,  Aberdeen  Proving  Ground,  MV, 
January  1965.  (AD  463926) 
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decelerating  effect  on  the  casing  is  less  than  would  be  predicted  by 
assuming  an  impulsive  start  at  the  time  the  shell  exits  the  gun. 

In  order  to  treat  in-bore  spin-up,  we  assume  that  Wedemeyor's  cc|ua- 
tion  based  on  impulsive  spin-up,  Eq.  (5),  can  be  used  in  cases  where  the 
spin  acceleration  is  very  large.  For  a typical  shell  the  spin  accelera- 
tion is  more  than  4,800  rev/s^,  making  it  almost  impulsive.  We  have  com- 
pared the  predicted  in-bore  angular  momentum  history  from  Method  I with 
results  from  a numerical  solution  of  the  Navier-Stokes  equations. 

In  our  work  the  in-bore  shell  motion  is  predicted  using  an  interior 
ballistics  trajectory  computer  program,  and  this  history  is  used  together 
with  the  gun  twist  to  determine  the  casing  spin  rate  boundary  condition. 
Figure  2 shows  an  example  of  a typical  in-bore  spin  history  for  an  XMb87 
shell  launched  at  charge  three  from  the  M109  howitzer  using  M.^  propel- 
lant. The  variable  t*  is  time  measured  from  the  beginning  of  projectile 
in-bore  motion.  The  spin  rate  increases  as  the  shell  accelerates  down 
the  tube,  exiting  19  ms  after  firing  with  a spin  rate  p^^  = 92, 7>  rev/s, 
corresponding  to  a muzzle  velocity  = 284.7  m/s. 

The  angular  momentimi  acquired  by  the  liquid  during  in-bore  spin-up 
depends  on  the  in-bore  spin  history,  the  cavity  dimensions  and  the  fluid 
properties.  In  Method  I the  in-bore  angular  momentum  history  is  calcu- 
lated using  Wedemeyer's  model  by  specifying  the  shell  spin  history  dur- 
ing in-bore  travel.  Figure  7>  shows  predictions  for  the  in-bore  .angular 
momentum  history  for  two  modified  XM687  shell.  Round  7670  contains  a 
highly  viscous  oil  (v  = 5x10”''  m^/s)  witl\  Re^^  = 5520;  the  laminar  endwall 

boundary- layer  assumption,  F;q.  f7)is  used.  This  calculation  predicts 
that  the  liquid  angular  momentum  is  18.6%  of  at  shell  exii.  Round 

7675  contains  water  (v  = Ixio"®  m^/s)  and  has  a much  higher  launch 
Reynolds  number  (Re^  = 1.7x10*^);  therefore,  the  turlnilent  endwall  l)ound- 

ary  layer  assumption,  F.q.  (8),  is  used.  The  water  payload  tlius  si^ins  up 
much  slower  than  the  oil  acquiring  only  2.6%  of  1-^  at  muzzle  exit. 

In  Method  III  a finite-difference  procedure  is  used  to  solve  the 
Navier-Stokes  equations.  The  flow  is  assumed  axisymmetric,  witli  u,  v, 
w and  F functions  of  r,  z and  t.  This  formulation  does  not  require  an 
endwall  compatibility  condition,  such  as  Eq.  (7).  T!\e  entire  flow  is 
calculated  without  separating  the  problem  into  boundary  layer  and  core 
regions.  Analytical  transformations  of  r and  z permit  grid  points  to  l)c 
densely  grouped  near  the  endwalls  and  sidewall  to  aid  in  resolving  the 
boundary  layers.  A numerical  solution  is  obtained  for  equations  des- 
cribing vorticity  and  circulation  using  a modification  of  the  predictor 
corrector  multiple  iteration  method  of  Rubin  and  The  stream 


6.  S.  G.  Rubin  and  T.  C.  Lin,  "A  Numevioal  Method  for  Threc-Dlmeneumal 
Viscous  Flow:  Application  to  the  Hijpevsoyiia  Leading  Fdge,  " J . Comp. 

Phys, , Vol.  9,  1972,  pp.  3,‘^9-3G4. 
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function  equation  is  solved  by  successive  line  relaxation.  A separate 
report^  will  give  details  of  the  Navier-Stokes  procedure  and  compare 
results  from  it  with  those  from  Wedemeyer's  model. 

Results  from  the  Navier-Stokes  equations  for  the  in-bore  angular 
momentum  history  are  shown  by  the  dashed  curve  in  Figure  3 for  Round 
7670.  This  calculation  used  a 21x41,  r-z  grid  and  took  22  time  steps  to 
reach  muzzle  spin  conditions,  requiring  5.2  minutes  of  CPU  time  on  the 
UNIVAC  1108.  Longer  computing  times  would  be  required  for  larger 
Reynolds  numbers.  Method  III  predicts  an  angular  momentum  level  that  is 
2.5%  larger  at  muzzle  exit  than  that  predicted  by  Method  I.  The  dis- 
crepancy between  these  two  calculations  remains  almost  constant  after 
about  6 ms.  A Navier-Stokes  calculation  has  not  been  made  for  Round 
7675  because  of  the  large  Reynolds  number  and  the  expected  turbulent 
flow  in  the  endwall  boundary  layer.  On  the  basis  of  the  comparisons  for 
Round  7670,  we  feel  justified  in  using  Method  I for  the  in-bore  angular 
momentum  history  of  the  liquid. 


IV.  COMPARISON  WITH  MEASURED  SPIN  DECAY 

Spin  decay  predictions  have  been  made  for  XM687  shell  launched  from 
the  M109  howitzer  and  compared  with  measurements  obtained  using  a solar 
aspect  sensor  and  telemetry  technique^.  The  motion  of  the  projectile  is 
determined  by  the  use  of  photovoltaic  cells  which  sense  the  orientation 
of  the  shell  relative  to  the  sun.  The  measuring  system,  called  a yaw- 
sonde,  is  carried  on-board  the  projectile  and  data  are  transmitted  to  a 
ground  station  through  an  FM/FM  telemetry  link.  Data  are  not  acquired 
during  the  in-bore  spin-up  process,  but  are  first  received  shortly  after 
the  shell  exits  from  the  gun  tube.  The  data  reduction  procedure  yields 
both  the  yawing  motion  of  the  projectile  and  the  spin  history  over  the 
whole  trajectory. 

The  standard  XM687  contains  two  liquid-filled  cylindrical  canisters 
separated  by  discs  which  rupture  on  launch,  producing  a single  cylindri- 
cal cavity.  Data  which  can  be  compared  directly  with  the  present  theory 
were  obtained  by  Mark®  at  Wallops  Island,  VA,  in  May  1975  using  a non- 
standard cavity  in  the  XM687  shell.  The  two-canister  configuration  was 
replaced  with  a single  cylindrical  cavity  having  a height  of  0.474  m and 
a diameter  of  0.107  m.  These  dimensions  are  approximately  the  dimensions 
of  the  cylindrical  cavity  in  the  standard  XM687  after  rupture  of  the 
discs.  Shell  filled  with  oil  and  with  water  were  fired  to  give  two 
quite  different  launch  Reynolds  numbers,  approximately  3.3x10®  and 


7.  C.  V.  Kitahens,  Jr.,  "Nauier-Stokea  Solutions  for  Spin-Up  from  Rest 
in  a Cylindriaal  Container,"  BRL  Report  (in  preparation) . 

9.  A.  Mark,  "Meaeurement  of  Angular  Momentum  Transfer  in  Liquid- Filled 
Shell,"  BRL  Report  (in  preparation). 
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1.7x10^,  respectively.  The  standard  XM687  is  approximately  90^o  filled, 
and  in  these  special  firings  fill  ratios  of  both  90  and  100%  were  used. 
Table  I provides  a summary  of  the  firing  data  for  the  five  rounds  for 
which  comparisons  are  presented  between  the  yawsonde  measurements  and 
the  predicted  spin  decay.  Round  10G2,  containing  a dual  canister,  was 
fired  at  Nicolet,  Canada,  during  the  1974-75  winter  tests^  under  dif- 
ferent meteorological  conditions  than  the  other  rounds.  For  all  calcu- 
lations discussed  in  this  report  the  spin  damping  coefficient  was  des- 
cribed by 

= 0.00860  M - 0.0200  for  0.53  < M < 0.84  , (19) 

P 

where  M is  the  instantaneous  free  flight  Mach  number.  liq.  (19)  was 
determined  from  firings  of  XM687  containing  solid  filler  as  discussed  in 
Reference  8.  The  air  density  and  temperature  needed  to  define  the  spin 
damping  function,  liq.  (13),  were  determined  from  meteorological  data  and 
the  projectile  velocity  was  deduced  from  radar  measurements.  The  spin 
damping  functions  for  two  representative  ca«es  are  shown  in  Figure  4. 
Round  10G2  was  fired  at  05®  quadrant  elevation  with  high  air  density. 
Round  7675  was  fired  at  30°  quadrant  elevation  with  almost  standard 
density.  All  rounds  exhibit  the  concave- doimward  shape  shown  in  Figure 
4.  The  spin  damping  functions  for  Rounds  7670,  7676  and  7677  arc  almost 
identical  to  that  for  Round  7675. 

A.  Fully-Filled  Shell 

The  predicted  spin-decay  for  round  7670  is  compared  with  the  yaw- 
sonde  measurement  in  Figure  5.  This  round  was  filled  with  oil 
(v  = 5x10“^*  m^/s) ; the  firing  data  for  this  and  subsequent  rounds  arc 
shown  in  Table  I.  The  total  time  of  flight  was  28s;  only  the  first  ten 
seconds  are  shown  in  Figure  5 because  the  spin-up  process  is  completed 
by  this  time.  The  spin-decay  history  deduced  from  the  yawsonde  measure- 
ment is  shown  by  the  solid  curve.  Data  acquisition  began  0.1s  after 
launch  and  the  measured  data  were  smoothly  extrapolated  hack  to  a calcu- 
lated value  of  spin  at  the  muzzle.  For  each  round  the  spin  at  the 
muzzle  was  calculated  using  the  known  gun  tube  twist  and  the  independ- 
ently measured  muzzle  velocity. 

In  the  theoretical  results  shown  in  F'igure  5 (note  the  broken 
ordinate  scale)  the  endwall  boundary- layers  were  assumed  to  be  laminar. 
Assuming  an  impulsive  start,  calculations  were  made  using  Method  I 
(shown  by  filled  circles)  and  Method  II  (shown  by  filled  triangles). 

The  results  from  these  two  calculations  are  almost  identical.  The 


9.  V.  Qokay  and  J.  H.  Whiteside,  "Flight  Behavior'  of  265mm  (XM687  Mod  I 
and  XM687  Mod  II)  and  8-Inah  (XM736  Mod  I)  Binary  Shell  at  Niaolet, 
Canada,  During  the  Winter  of  1974-1976, " BRL  Memorandum  Report  2608, 
Aberdeen  Proving  Ground,  MD,  Maroh  1976.  (AD  B010566L) 
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Table  I.  Firing  Data  for  XM687  Shell 


calculated  spin  decay  rate,  assuming  an  impulsive  start,  is  larger  than 
the  observed  decay  rate.  This  is  due  to  the  neglect  of  the  in-bore 
spin-up  process.  At  t = 10s  the  difference  is  0.4%. 

When  in-bore  effects  are  included  in  Method  1,  the  results  (shown 
by  open  circles  in  Figure  5)  agree  better  with  the  measurement.  The 
in-bore  spin-up  process  described  in  Section  III.  C was  used  for  the 
19ms  prior  to  t = 0.  Time  zero  is  defined  to  be  the  time  when  the  pro- 
jectile clears  the  gun  tube.  The  predicted  spin-up  process  was  com- 
pleted at  t = 1.4s;  for  t > 1.4s  the  shell  is  in  the  late  free-flight 
phase,  shown  in  Figure  Ic,  wherein  the  direction  of  the  liquid  moment 
reverses.  At  this  time  the  predicted  spin  rate,  including  in-bore 
effects,  was  0.13%  larger  than  the  measured  value.  Calculations  includ- 
ing in-bore  effects  have  not  been  n>ade  for  Method  IT.  The  predicted  and 
measured  spin-decay  curves  cross  each  other  at  18s  and  the  predicted 
spin  rate  is  0.43%  lower  than  the  measured  value  when  the  yawsonde  data 
end  at  26.6s.  This  cross-over  can  probably  be  traced  to  inaccuracies  in 
the  measured  which  causes  the  spin  damping  term  f(t)  to  be  too 
large.  p 

When  the  launch  Reynolds  number  is  greater  than  3x10^  the  endwall 
boundary  layers  arc  expected  to  be  at  least  partially  turbulent.  The 
measured  spin  decay  for  round  7675,  filled  with  water,  is  compared  with 
results  from  the  present  theory  with  turbulent  Fikman  layers  in  Figure  6. 
Methods  I and  II  both  predict  a more  rapid  spin  decay  than  is  observed 
experimentally;  Method  I gives  the  better  agreement.  At  this  large 
Reynolds  number  in-bore  effects  are  not  significant  and  both  calcula- 
tions with  Method  1 (filled  and  open  circles)  give  similar  results.  The 
predicted  time  of  reversal  of  the  liquid  moment  is  24.5s.  At  this  time 
the  predicted  spin  rate  is  0.79%  lower  than  the  measurement,  as  extrapo- 
lated from  the  available  yawsonde  data  which  ends  at  23.9s. 

B.  Partially-Filled  Shell 

At  the  present  time  a spin-up  model  for  the  partial ly- fi 1 led  case 
does  not  exist,  and  the  treatment  of  many  cases  of  interest  is  hampered 
by  the  lack  of  such  a theory,  since  most  liquid-filled  shell  are  only 
partially-filled.  It  would  be  useful  to  learn  if  the  methods  developed 
here  for  the  filled  cylinder  could  be  used  with  any  confidence  for  pre- 
dicting spin  decay  in  partially-filled  cylinders.  We  might  expect  this 
if  the  fill  ratio,  B,  is  not  much  less  than  90%. 

The  presence  of  the  free  surface  in  the  partially-filled  case 
introduces  many  complexities  to  the  spin-up  flow  during  in-bore  acceler- 
ation and  launch.  '’Tiis  complex  motion  is  discussed  in  Reference  10. 


1 0,  Snginaering  Design  Handbook^  Liquid- Filled  Projectile  Deoigrit  AMC 
Pamphlet  No,  706-165^  U.S.  Amy  Materiel  Command,  Washington,  DC, 
April  1969,  p.  8-4.  (AD  853719) 
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Briefly,  the  liquid  moves  rearward  in  the  cavity  under  the  action  of  in- 
bore accelerating  forces  and  then  forward  due  to  air  drag  as  the  shell 
enters  free  flight.  The  process  whereby  the  free  surface  forms  near  the 
rear  endwall  of  the  cavity  as  the  shell  emerges  from  the  gun  tube  may 
vary  considerably,  depending  on  the  actual  fill  ratio  of  the  cavity. 

This  liquid  motion  in-bore  and  during  the  early  part  of  the  flight 
possibly  aids  in  mixing  the  rotating  fluid  near  the  cavity  sidewall  with 
the  non-rotating  fluid  away  from  the  wall.  The  importance  of  this 
effect  on  the  spin  decay  history  of  partially-filled  shell  is  not  pres- 
ently known. 

The  Wedemeyer  spin-up  theory  has  been  applied  to  partially-filled 
XM687  shell,  treating  the  shell  as  though  it  were  fully-filled.  The 
spin  history  of  two  XM687  (again  having  a single  cylindrical  cavity) 
with  B = 0.90  are  compared  with  predictions  using  this  theory.  Compari- 
sons for  Round  7677  (Table  I)  are  shown  in  Figure  7.  The  spin  history 
for  this  laminar,  90%-filled  case  is  almost  identical  to  that  observed 
for  the  fully-filled  case  shown  in  Figure  5.  Methods  1 and  II  both  give 
reasonable  agreement  with  the  measured  spin  decay.  Method  I,  including 
in-bore  effects,  is  slightly  more  accurate;  at  10s  it  predicts  a spin 
rate  that  is  0.18%  larger  than  the  yawsonde  measurement,  whereas  Method 
II  predicts  a value  that  is  0.29%  smaller. 

Similar  comparisons  for  a much  higher  Reynolds  number  are  shown  in 
Figure  8 for  round  7676  which  is  90%-filled  with  water.  It  is  assumed 
that  the  endwall  boundary  layers  are  turbulent  in  this  case.  The  spin 
history  is  qualitatively  different  than  that  observed  for  the  100%- 
filled  case  in  Figure  6.  The  most  noticeable  difference  occurs  in  the 
first  few  seconds  of  the  flight;  the  90%-filled  case  loses  spin  much 
faster  than  the  100%-filled  case.  In  the  first  half-second  of  flight 
the  90%-filled  case  loses  2.0%  of  its  muzzle  spin  rate,  whereas  the 
100%-filled  case  only  loses  1,2%.  This  is  surprising  because  it  seems 
to  indicate  that  the  spin-up  process  in  the  90%-filled  case  with  a 
turbulent  endwall  boundary  layer  is  faster  than  in  the  100%-filled  case. 
This  may  possibly  be  due  to  the  mixing  caused  by  the  liquid  shifting 
from  the  rear  to  the  front  of  the  cavity  as  the  round  emerges  from  the 
gun  tube.  Definite  conclusions  about  this  cannot  be  reached  because  B 
is  not  the  only  parameter  that  varies  in  these  two  firings. 

Figure  8 shows  that  for  3 = 1.00  and  turbulent  Ekman  layers,  the 
theory  predicts  a more  gradual  spin  decay  in  the  first  few  seconds  of 
the  flight  than  is  observed  for  the  B ® 0.90  round.  Method  I gives  a 
more  accurate  prediction  of  the  shell  spin  decay  than  Method  II.  At  10s 
Method  II  predicts  a spin  rate  which  is  in  error  by  approximately  1% 
while  Method  I essentially  matches  the  experimental  measurement.  Since 
this  shell  with  B = 0.90  doesn't  appear  to  have  the  same  spin  decay 
history  as  the  100%-filled  shell  one  should  be  cautious  in  attempting  to 
apply  Wedemeyer 's  fully-filled  spin-up  theory  to  this  case. 


ir;ninu  :;.ig^j4n3i4, . 
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V.  DISCUSSION  OF  FLUID  SPIN-UP  TIMi; 


Solid-body  rotation  is  never  completely  achieved  in  a liquid-filled 
projectile  because  p is  constantly  changing.  This  is  illustrated  in 
Figure  9 which  shows  liquid  azimuthal  velocity  profiles  predicted  with 
Method  I at  three  instants  in  the  flight  of  round  7675.  The  azimuthal 
velocity  normalized  by  ap^  at  the  sidewall,  r/a  = 1.0,  decreases  with 

time.  A typical  velocity  profile  for  solid-body  rotation  is  shown  for 
t = 27.5s  by  the  dashed  straight  line  v^  = rp.  This  state  is  not 

achieved  at  any  time  during  this  flight;  at  t = 27.5s  v^  < rp  for 

0 < r/a  < 0.57  and  v^  > rp  for  0.57  < r/a  < 1.  The  fluid  in  the  latter 

range  of  r is  "overspun"  relative  to  the  projectile  at  this  time;  i.e., 
the  local  fluid  angular  velocity,  = v^/r,  is  greater  than  p.  Shortly 

after  the  liquid  near  the  sidewall  becomes  overspun  reverses  direc- 
tion, opposing  The  spin-up  process  is  then  in  the  late  free- 

f light  phase  shown  in  Figure  Ic. 

At  a given  time,  solid-body  rotation  is  more  closely  approached  for 
small  Re^,  because  diffusion  effects  are  greater.  Tins  is  illustrated 

in  Figure  10  which  shows  predicted  for  round  7670.  The  approach 

to  solid-body  rotation  for  roi|nd  7670  is  much  more  rapid  than  for  round 
7675,  primarily  due  to  the  mudh  lower  Reynolds  number.  Also,  the  liquid 
in  round  7670  approaches  .soli|-body  rotation  much  more  closely.  As  an 
example,  at  t = 27.5s  the  largest  difference  between  the  predicted  azi- 
muthal velocity  and  solid-body  rotation  is  less  than  0..5()%;  for  round 
7675  differences  of  4o  for  r/a  - 0.90  and  -IS'l,  for  r/a  = 0.20  are 
obtained. 

It  would  be  convenient  to  have  a unique  definition  of  spin-up  time, 
but  this  is  not  possible  because  of  the  asymptotic  approach  to  solid- 
body  rotation.  Also,  depending  on  the  application,  different  measures 
of  the  approach  to  solid-body  rotation  are  appropriate.  A comparable 
situation  exists  in  defining  the  thickness  of  a boundary  layer.  Since 
this  report  is  concerned  with  spin  decay,  we  adopt  a definition  which 
uses  the  calculations  of  Method  I.  A separate  report  will  discuss  other 
definitions. 


The  spin-up  time,  t^,  is  the  time,  measured  from  t = 0,  at  which 

the  fluid  angular  momentum  reaches  99%  of  the  rigid-body  angular  momen- 
tum, when  the  latter  is  based  on  the  instantaneous  shell  spin  rate  p. 
Figure  11  shows  angular  momentum,  L,  for  rounds  7670  and  7675  as  pre- 
dicted by  Method  I,  including  in-bore  effects.  L is  normalized  by  L^, 

the  angular  momentum  the  fluid  would  have  if  it  were  rotating  as  a 
rigid-body  at  the  instantaneous  shell  spin  rate.  As  t increases 
L/L^  ->•  1;  it  can  exceed  unity  late  in  the  flight  when  the  fluid  is  in 


the  overspun  state.  This  can  occur  even  if  v^  < rp  for  fluid  at  small  r 

because  the  overspun  fluid  at  large  r gives  the  dominant  contribution  to 
the  integral  for  L,  Eq.  (11).  Late  in  the  flight  L/L^  > 1 for  both 

rounds  in  Figure  11,  but  it  is  only  obvious  for  round  7675  because  of 
the  small  scale  of  the  plot.  Normalizing  L with  the  constant  would 

be  convenient,  but  then  L/L^  is  non-monotonic  and  has  a maximum  that  is 

always  less  than  one.  According  to  our  definition  of  spin-up  time,  we 

obtain  t = 0.9s  for  round  7670  and  t » 19.0s  for  round  7675.  The 
s s 

large  difference  in  values  of  t^  results  from  the  factor  of  10^  differ- 
ence in  Reynolds  number. 

It  should  be  pointed  out  that  the  angular  momentum  predictions 
shown  in  Figure  11  do  not  agree  with  the  results  shown  in  Figure  13  of 
Reference  2,  wherein  the  instantaneous  liquid  angular  momentum  is  cal- 
culated by  fitting  measured  yawsonde  spin  rate  data  and  numerically 
integrating  the  projectile  roll  equation.  For  a typical  XM687,  round 
El-5977  with  Re^  = 1.7x10^,  it  is  concluded  in  Reference  2 that  only  85% 

of  the  instantaneous  rigid-body  angular  momentum  is  achieved  in  the  30s 
time  of  flight  The  present  calculations  predict  that  85%  of  is 

achieved  in  approximately  0.7s  for  round  7675  with  Re^  = 1.7xio^.  The 

large  difference  between  these  results  for  spin-up  time  appears  to  be 
caused  by  the  inaccurate  value  for  C»  used  in  Reference  2 and  the 

P 

neglect  of  in-bore  effects®.  More  recent  calculations  of  angular  momen- 
tum transfer®  show  better  agreement  with  the  present  results. 


VI.  CONCLUSIONS 


This  report  has  described  two  methods  based  on  Wedeineyer's  spin-up 
theory^  for  predicting  the  spin-decay  of  liquid-filled  projectiles 
throughout  their  flight.  Shell  spin  histories  are  predicted  by  coupling 
the  solution  of  the  axial  spin  decay  equation  for  the  projectile  to  the 
solution  for  the  liquid  spin-up  process. 


Method  II,  based  on  an  ordinary  differential  equation  for  liquid 
angular  momentum,  is  strictly  applicable  only  during  the  spin-up  pro- 
cess. It  does  not  apply  late  in  the  flight  when  the  liquid  is  "over- 
spun"  relative  to  the  shell  spin  rate;  i.e.,  the  shenr  stress  acting  on 
the  casing  reverses  direction  and  opposes  further  shell  spin  decay. 

This  effect  is  accounted  for  in  Method  I,  based  on  a numerical  solution 
of  Wedemeyer's  spin-up  equation,  Eq.  (5).  Here  the  diffusion  term  is 
retained  in  Wedemeyer's  equation  and  finite-difference  solutions  are 
obtained  for  the  spin-up  process.  Method  I gives  better  agreement  with 
experimental  measurements  for  shell  spin  decay  than  Method  II. 
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The  liquid  angular  momontum  in-bore  has  been  calculated  using 
numerical  solutions  of  the  Navier-Stokes  equations  and  Wcdemeycr's  eciua- 
tion.  The  results  show  that  in  some  XM687  firings,  with  high  viscosity 
oil,  approximately  IS"®  of  the  liquid  rigid-body  angular  momentum  is 
attained  in-bore  in  less  thar  one  revolution  of  the  shell.  The  inclu- 
sion of  in-bore  spin-up  in  the  spin  decay  calculations  with  Method  1 
improved  the  agreement  with  the  experimental  measurements,  significantly 

more  so  for  small  Re  than  large  Re  . 

o * o 

Calculations  for  XM687  shell  spin  decay  were  compared  with  yawsonde 
measurements  for  100%  and  90%-filled  cases.  The  100%-filled  cases  show 
excellent  agreement  between  the  theory  and  the  measurement  throughout 
the  entire  flight,  the  errors  being  less  than  0.8%  of  the  measured  spin 
rate.  Spin  decay  for  90%-fillcd  cases  are  not  quite  as  accurately 
described  by  the  present  theory,  but  the  errors  are  no  larger  than  about 
1%.  The  largest  errors  occur  in  the  90%-filled  case  at  high  Reynolds 
number  where  the  boundary  layers  are  expected  to  be  turbulent  over  part 
of  the  endwall. 

The  spin  records  discussed  in  this  report  all  exhibit  smooth  spin 
decay  throughout  the  flight.  There  are  other  types  of  spiii  records  that 
cannot  be  predicted  by  the  present  theory.  When  the  XM687  slicll  is 
launched  with  large  yaw  the  spin  history  can  be  quite  different  from 
that  shown  in  figures  5-8.  In  these  cases  the  spin  history  is  char- 
acterized by  a much  more  rapid  initial  spin  decay,  which  is  almost 
linear  for  approximately  one  to  two  seconds,  followed  by  a sharp  change 
to  a smaller  rate  of  decay.  This  spin  record  with  a "corner"  is  not 
predicted  by  the  methods  described  in  this  report.  A method  for  treat- 
ing the  spin  decay  in  such  cases  is  discussed  in  Reference  3. 

Calculations  with  Method  I have  been  used  to  determine  the  spin-uip 
time,  t^;  defined  as  the  time  when  99%  of  the  instantaneous  rigid-body 

angular  momentum  is  achieved,  for  the  XM687  considered  here,  t_,  varies 

from  0.9s  for  Re  = 3320  to  about  19s  for  Re  ~ 1x10^’.  These  predictions 

for  t^  are  based  on  Wedemeyer's  spin-up  model  winch  is  valid  for  small 

yaw.  When  the  projectile  yaw  is  large,  the  toroidal  vortex  model  dis- 
cussed in  Reference  3 may  be  more  appropriate  for  determining  a measure 
of  the  spin-up  time. 
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a.  IN-BORE  b.  EARLY  FREE  FLIGHT  C.  LATE  FREE  FLIGHT 


Figure  1.  Moments  Acting  on  a Spinning  Liquid-Filled  Projectile; 
A Free-Body  Diagram  of  Casing  During  Launch  and 
Free  Flight 


Figure  2.  Normalized  In-Bore  Spin  History  for  a Typical  XM687 
Shell  (Round  7670)  Launched  from  the  Ml 09  Howitzer 
at  Charge  3,  p = 92.3  rev/s 
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NAVIER-STOKES  SOLUTION  . 

WEDEMEYER  MODEL  (METHOD  I)  // 

// 

// 

// 

// 

// 

// 

// 

// 

//Ground  7670 

/ Reo=  3320 

// 

// 

//  ROUND  7675 
//  Re«=1.7xl0i 


Figure  3.  Normalized  In-Bore  Liquid  Angular  Momentum  History  for  Two 
XM687  Shell  with  Different  Launch  Reynolds  Numbers 


■ROUND  IOG2 


.ROUND  7675 


Figure  4.  Spin  Damping  Function  f(t)  for  Two 
Representative  XM687  Shell 


t(s) 


Figure  5.  Shell  Spin  Decay  for  100%-Filled  Round  7670  with 
R6q  = 3320  and  Laminar  Endwall  Boundary  Layer 
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Figure  7.  Shell  Spin  Decay  for  90%-Filled  Round  7677  with 
Re^  = 3421  and  Laminar  Endwall  Boundary  Layer 
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SHELL  SPIN  RATE  AT  MUZZLE 
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Figure  9.  Azimuthal  Velocity  Profiles  for  Round  7675  with 

Re^  = 1.7x10®  and  Turbulent  Endwall  Boundary  Layer 


Figure  10.  Azimuthal  Velocity  Profiles  for  Round  7670  with 
Re„  = 3320  and  Laminar  Endwall  Boundary  Layer 
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LIST  OF  SYMBOLS 


a radius  of  shell  cylindrical  cavity  (=  0.0535  m for  non-standard 

XM687  cavity) 

c half-height  of  cylindrical  cavity  (=  0.237  m for  non-standard 

XM687  cavity) 

f(t)  defined  by  Eq.  (13)  (kg-m^/s) 
t projectile  diameter  (=  0.1524  m for  XM687  shell) 

p instantaneous  spin  rate  of  shell  (s“^) 

p^  shell  spin  rate  at  muzzle  (s“^) 

r radial  coordinate  (m) 

t time  measured  from  instant  of  projectile  release  from  gun 

tube  (s) 

t*  time  measured  from  the  beginning  of  projectile  in-bore 
motion  (s) 

t time  at  which  the  fluid  angular  momentum  reaches  99-5  of  the 

rigid-body  angular  momentum  based  on  p (s) 

« 

t 

u,  radial  velocity  component  in  core  fl(?w  in  Wcdcmeycr  spin-up 
model  (m/s)  / 

/ 

V azimuthal  velocity  component  in  liquid  (m/s) 

/ 

/ 

v azimuthal  velocity  component  in  core  flow  in  Wedemeyer  spin-up 

model  (m/s) 

w,  axial  velocity  component  in  core  flow  in  Wedemeyer  spin-up 
model  (m/s) 

z axial  coordinate  measured  from  cavity  rear  ondwall  (m) 

projectile  roll  damping  coefficient  (nondimensiunal) 

P 

axial  moment  of  inertia  of  projectile  casing  (kg-m'^) 

L liquid  angular  momentum  (kg-m^/s) 

L^  rigid-body  angular  momentum  of  liquid  at  spin  rate  p^  (kg-m^/s) 

L^  rigid-body  angular  momentum  of  liquid  at  spin  rate  p (kg-m^/s) 


I 
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LIST  OF  SYMBOLS  (continued) 


’IB 


M 


Aero 


M 


Liq 


P 

c 

Re 


Re 

o 


in-bore  moment  acting  on  shell  due  to  gun  tube  twist  (N-m) 

aerodynamic  moment  acting  on  shell  due  to  air  shear  (N-m) 

liquid  moment  acting  on  shell  due  to  shear  stress  at  cavity 
walls  (N-m) 

pressure  (N/m^) 

Reynolds  number  of  shell  defined  by  Eq.  (9)  (nondimensional) 
launch  Reynolds  number  of  shell  defined  by  Eq.  (14)  (nondimensional) 


S maximum  cross-sectional  area  of  the  projectile  (=  0.01824  m^  for 
XM687  shell) 

muzzle-velocity  of  projectile  (m/s) 

V instantaneous  projectile  velocity  (m/s) 

3 fill  ratio  of  cylinder;  fluid  volume  divided  by  total  cavity 
volume  (nondimensional) 

e tolerance  on  angular  momentum  convergence  from  iteration  to 

iteration 

0 angular  coordinate 

V kinematic  viscosity  of  liquid  (m^/s) 

p liquid  payload  density  (kg/m^) 

p air  density  (kg/m^) 

£i 

0)  local  liquid  angular  velocity  (=  v/r)  (s"^) 

0)  local  liquid  angular  velocity  in  core  (=  v /r)  (s"^) 

c c 

Superscript 

1 iteration  level  i 


Subscript 

n time  level  n 
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